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1 LI 

A METHpD AND A SYSTEM FOR CARDIAC ANALYSIS 

The invention relates to a method according to the proamble of the 
appended independent claim 1 for cardiac analysis, which method 
comprises steps for acquiring an ECG-signal, detecting at least one 
wave of the ECG-signal and calculating parameter values of said wave. 
The Invention also relates to a cardiac analysis system according to the 
preamble of the appended independent claim 18 for Implementing the 
aforementioned method as well as to a computer program product 
according to the preamble of the appended Independent claim 30. 

: Background of the invention 

An electrocardiogram (ECG) is a record of electrical events occurring 
15 within the heart of a patient and it is used for examination of heart 
diseases. The ECG plots the sequence of changes in electrical 
potential differences recorded between different regions of the patient's 
hody surface. The ECG can be measured using various lead systems 
known as such. Commonly the ECG is obtained by using a standard 
20 12-lead arrangement, but it can be obtained by using other lead 
systems, for example a Frank system, which uses three leads. These 
ECG-measuremenls provide a one-dimensional electrocardiographic 
signal, but ft is also possible to present ECG as a three-dimensional 
signal, a so-called vectorcardiography. 

25 

The basic form of the ECG-signal can be seen in Figure 1. The ECG 
breaks down each heartbeat into a series of electrical waves. Three ot 
the waves (P. QRS, T in the figure), a P-wave. a QRS-complex and a 
T-wave, are associated with the hearts contractions. The first 

30 deflection of the ECG is called the P-wave. It reflects the sequential 
activation (depolarization) of the heart's upper chambers, the right and 
the left atria. After the P-wave the ECG returns to its baseline, which is 
a substantially straight line on an ECG-paper where there are no 
positive or negative changes of electricity to create deflections. The 

35 next Is a complex called QRS-complex, which reflects the sequential 
activation of the heart's lower chambers, the right and the loft ventricle. 



5 

V 

10 



VAST OTTO 10-04-2003 1B:07 MISTA- 03 2886262 



KENELLEPATREK Asiakaspalvel SIVU 002 



2 



After the QRS-complex, tho ECG returns to, or very nearly to, its 
baseline which is nailed a ST-segment, where il remains until' the 
appearance of a T-wave. Tho T-wavo reflects ventricular repolarization. 

5 When discussing the ECG, a term "interval" refers to the length of the 
wave along with an Isoelectric line following it. The interval can be 
named by using the letters of both waves on either side (e.g. PQ) A 
space between dash lines 1 and 3 (Figure 1) represents the PQ- 
intcrval. A term "segment" refers to the baseline between the end of 
10 one wave and the beginning of the next wave. The space between 
dash lines 2 and 3 represents the PQ-segment, as well as the space 
between dash lines 4 and 5 represent the ST-seymenl. 

The above-mentioned vectorcardiogram (VCG) Is a form of the 
15 electrocardiography that represents movements of vectors of the 
myocardial activation with direction and magnitude as loops in a three- 
dimensional space. Vectorcardiography can be considered as a 
method or recording the direction and magnitude of electrical forces of 
the heart by means of a continuous series of vectors that form a 
20 curving line around a center. A spatial orientation and magnitude of the 
heart vector is projected onto three orthogonal planes as frontal, 
horizontal and sagittal planes. The tip of the vector with coordinates (X, 
Y, Z) should trace out a loop in space, called a vector loop. The vector 
loop starts from the zero-point, which corresponds to the isoelectric line 
25, or the baseline on a usual scalar ECG. The loop consists of, similar to 
the ECG-signal, a P-loop, a QRS-Ioop and a T-loop. When heart cycle 
is completed, the QRS-loop returns, in a normal situation, to ?ero-pnint 
and closes the loop. On the other hand, if the patient has a ST-change, 
the QRS-complex will not end to the baseline, but over or under 
30 depending on whether the patient has a ST-e!evatlon or depression. ' 

The general features of the QRS-complex and the ST-segment wHhin 
the FCG-signal have been exrensiveiy S Turiieri and their diagnostic 
significance is quite well known. On tho contrary, far less attention has 
35 been paid to the P-wave and therefore the clinical properties of Ihis 
waveform are still not fully understood. It is known that the P-wave 
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represents atrial depolarization and it can show the heart's rate and 
rhythm. It is known that the atrial abnormalties and the P-waves are 
related to some heart diseases, tor example to heart failure or to the 
susceptibility to atrial arrhythmias. Still, the specificity of these changes 
5 has been low. 

Developed ECG-systems observe the QRS-compiex and the ST- 
segment ideally, but have a lack in their operation of observing the P- 
wave. A reason for that can be found, in common opinion, where the P- 

10 wave and the atrias have been considered less meaningful compared 
to the QRS-wave, the ST-segment and the vontriclcs. Due to this, 
there are significantly less studies concerning the P-wave. In addition, 
the known studies of the P-wave use mostly the known 12-lead ECg! 
Suoh a 12-lead ECG-arrangement has a disadvantage in expression of 

15 the waves; it does not give a good three dimensional picture of the 
electrical waves, which neglects the small waves, for example the P- 
wave. Studies done of the P-wave utilize mostly the ECG recording 
system that expresses the heart situation In one moment, but is not 
capable of continuously showing the in-timo coming dynamic changes 

20 of the examined P-wave. 

Short dasi:rirrtin r i of |)m IhwiuHmm 

The present invention relates lu a cardiac analysis method and system. 

25 which takes into account said lacks in existing systems. More precisely,' 
Hie cardiac analysis melhod according lu the invention comprises the 
steps for acquiring the ECG-signal, detecting at least one wave of the 
ECG-signal and calculating the parameter values of said wave, which 
is the P-wave. The cardiac analysis method according to the invention 

30 is then focused to the dynamic changes of the P-wave. The cardiac 
analysis system according to the invention comprises means tor 
acquiring the ECG-signal, means for detecting at least one wave of the 
ECG-signal and means for calculating the parameter values ot said 
wave, which is the P-wave. The system is additionally adapted to focus 

35 thR analysis tn thR dynamic changes nf the P-wave. wherein said 
system also comprises means for comparing substantially every 
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detected P-wave to the reference P-wavc in defined time period. The 
computer program product according to the Invention comprises 
computer readable instructions for detecting at least one wave from the 
ECG-signal, which said wave Is a P-wave, whereupon Ihe computer 
program code comprises computer readable instructions for focusing to 
the dynamic changes of the P-wave, wherein said computer program 
code additionally comprises computer instructions configured to 
compare substantially every detected P-wave to a reference P-wave in 
defined time period. 



In accordance with tho invention, tho ECG-signal is processed in the 
form of vectorcardiogram that can be described by three orthogonal 
leads: X, Y, Z. P-wavc is at first recognized by a template and then the 
signals of the P-wave are averaged to form a smoolh wave form, which 
15 is used to determine the P-wavc specific parameter values 
continuously from the whole ECG-recording. The cardiac analysts 
method according to the invention is then addressed to the dynamic 
changes of said P-wave. That means that it addresses lo ihe changes 
of P-wave configuration that happen smoothly in the follow-up time. 
The method comprises means for comparing substantially every 
detected P-wave to the reference P-wave in a defined time period. The 
resulls from each parameter are presented as points on a trend graph. 
The oardiao analysis method aocording the invention also provides 
means lo notice and pick up the atrial exlrasystoles, to separate them 
into different sub-classes based on their morphology and then to 
average and analyze them in the same way as the P-wave. 

The analysis of the P-wave according to the invention, brings new tools 
to the analysis of atrial activation; in consequence the P-wave is 
considered to be more meaningful than before. With the method 
according to the invention, there is a first-time chance to see the 
dynamic changes of P-wave in a patient who has an acute myocardial 
infarction that affects also the atrial tissue. The atrial infarction is the 
entity that is supposed to affect the patient outcome by causing atrial 
arrhyimlas and rupiures of the atrial free wall In special situations Still 
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the dynamic changes in the P-wave configuration and the diagnostic 
criteries of the atrial infarction are unclear. 

The analysis of the P-wave according to the Invention provides also an 
5 effective way to monitor atrial arrhythmias, especially the paroxysmal 
atrial fibrillation. I he atrial fibrillation is the most common cardiac 
arrhythmia, the prevalence of which is increasing with the aging of the 
population. Because of its clinical importance and the lack of highly 
satisfactory management approaches, it is the subject of active clinical 

10 and research efforts. With the method according to the Invention, it is 
possible to observe the dynamic changes of P-wavo aftor the 
cardioversion when the normal sinus rhythm has been achieved and to 
try to find explanations for why some of the patients will maintain sinus 
rhythm after the cardioversion and some of them will not. Additionally. 

1 5 with the method according to the invention, it is possible to observe the 
dynamic changes of P-wave in situations where the risk of atrial 
fibrillation is high, for example after heart surgery. Because the method 
of the invention is based on dynamic changes of the P-wave, it also 
gives a possibility of observing the actions of treatment and medicine 

20 that is given to maintain sinus rhyrhm. 

The analysis of the P-wave according to the invention also provides an 
effective way to monitor the dynamic changos of P-wave in the patient 
who has an acute heart insufficiency. These P-wave abnormalities are 
25 discussed in more detail later in the description. It is obvious that with 
the method according to the invention, it is possible lo observe Ihe 
dynamic changes of P-wave also in other situations where the atrias 
are in the changing abnormal situation that affecte liie P-wave 
configuration as well. 

30 

The preferred embodiment of the invention is set forth in the drawings 
in the deLailed description which follows, and in the annexed claims! 
Further objects and advantages of the invention are also considered in 
the description. The Invention itself is defined with particularity in the 
35 claims. 
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Brief ctesr.rint|on nf th* drawings 

F,9Ure1 is a graphical representation illustrating the 

ECG-signal, 

H 9 ure 2 illustrates a flow chart of the steps of the method 

according to the invention, 

n 3 ure 3 illustrates a flow chart of the stops of the P-wave 

analysis according to the Invention, 

Figure 4a-e illustrates the P-wave in one-dimensional 

environment, 

15 Figure 5a-b illustrates the P wavo vector loop in two- 

dimensional environment, 



10 



20 



25 



Figure 6a-b Illustrates the P-wave vector loop in three- 

dimensional environment, 

Rguro 7 illustrates the P-wave in magnitude ECG, and 

Figurc8a-b illustrates the change area, the change vector 

angle and magnitude. 

Pstalled (iHscrirriifin of i\m t ir^f errH( j embogimgnj 

Figure 2 illustrates the steps of the method according to the invention 
as a flow chart. The present invention exploits the known three- 
30 dimensional electrical ECG-model called vectorcardiogram. ECG- 
signal acquisition (201) can be done in various ways. A basic method is 
to use eight standard ECG-surface electrodes, which are placed on the 
patient according to the Frank electrode system. I he electrodes are 
then used to form an CCG-vectorcardiogram by method known in prior 
art which FCG-vflotoroarriingram can he described by three orthogonal 
leads X, Y, Z. It is also possible to acquire the ECG signal by using a 
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standard twelve-lead-ECG -arrangement, which is stored and derived to 
form a vectorcardiogram and then further analyzed. According lo the 
invention, the ECG-signal is acquired preferably from existing 
vectorcardiogram data collected by another system (200). For 
5 example, the ECG-signal can preferably be acquired from the MIDA 
data storage unit. The MIDA Is a commercial monrtoriny system (sold 
by Philips Medical) made to analyze the ischemic changes of QRS- 
nomplex and ST-segment of EGG and it is widely used in hospitals. 
MIDA registers the eleotrical signals of the heart using said Frank 
10 electrode system and constructs the three-dimensional electrical model 
from thorn. The action of MIDA is discussed in more detail in US- 
5520191. It is obvious that the ECG-signal in the present invention can 
bo acquired by not only the mentioned but any known methods. 

15 After the acquisition, the ECG-signal is preproceseed (202) to minimize 
electrical artefacts. Preferably the raw data is filtered to remove the 
noise and to improve its signal-to-noise-ratio (SNR). The noise is 
filtered by some known signal processing method, which are not 
discussed more in this as ft can be considered obvious to the man 

20 skilled in signal processing. 

The preprocessed signal is according to the invention analyzed to 
detect the peaks of R-waves (203). A gap between two peaks is 
referred to by a term "beat". When the peaks are found, a time period 

25 between them is measured and the beat is stored. Beats are classified 
depending on their duration (time period). If the time between two 
consecutive R-waves changes suddenly under predetermined time, It 
can be assumed that the beat consists of an atrial extra systole (204) 
In that case the beat is stored to the atrial extra systole database for 

00 further analysis (205). The atrial extra systoles can be categorized into 
sub-classes according to their morphology. Afterwards, they are 
averaged and analyzed. The result of this examination will gi ve 
information of the changes of the atrial extra systoles as well as the 
number of the atrial extra systole sub groups. By picking the extra 

35 systoles from the ECG-signal. the resulting P-wave data will become 
as homogenous as possible, which is preferable to the analysis. There 
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is also a possibility that the change results from the ventricular extra 
systole when there is no atrial extra systole preceding the Uevianl- 
looking QRS-complex. In this case, the changed boat is also picked to 
the second database and later removed. All the uLher beats (the 
5 unchanged, time period within the predetermined time limit) are stored 
to a P-wave database. The beats In said both databases are stored in 
the X, Y, Z -leads (206) and used separately for the P-wave analysis 
according to the invention. 

10 P-wave analysis: 

Reference is now marie to Figure 3, in which the examination of the P- 
wave is described in more detail. At first the P-wave signal is 
processed by baseline correction (307) Lo minimize the baseline drift of 
15 the signal. This way the signal quality will improve. There are several 
methods for implementing the baseline correction, for example linear 
interpolation, quadratic or cubic corrections, polynomial fitting and high 
pass frequency filtering. Because said methods are known as such, 
they are not explained further. 

20 

The P-wave is normally viewed as a smooth, small and curved 
deflection. It is usually a positive deflection, although it may be 
negative as well. The duration of the P-wave is normally below 0.12 
seconds and the amplitude is normally less than 0.25 mV. In the frontal 
25 plane the P-wave axis angle ranges from 0 to /b degrees. Also 
notched P-waves may be seen and the normal P-wave may often 
exhibit two components as an M-shaped complex. Because of its 
small, ill-defined and variable shape, the detection of the P-wave is 
difficult. 

30 

According to the invention, the P-wave detection (308) is preferably 
done by a template method, which is a type of two-dimensional cross- 
correlation method with both time and amplitude as variables. The 
template method attempts to detect the P-wavc by using the 
35 covarlance between the template signal and the wave form signal 
stored boat For substantially each beat, the template is tried to be 
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matched with the actual signal and the P-wave is detected when the 
covariance exceeds a specifier! threshold. The P-wave is valid if its 
duration is between the minimum and the maximum P-wave durations. 
As said, the length of the P-wave is generally about 120 ms. Time 
between the beginning of the P-wavc to tho beginning of the Q-wave 
(called PQ-tlme) is generally below 200 ms. These values are 
preferably used to form a template time window, but the time window 
can be predefined with some other values or with entering the values 
occasionally. 



The P-wave detection is preferably carried out only in one coordinate, 
because frequently one coordinate is less noisy than Ihe others and 
therefore more suitable for analysis. In the detection, the user is at first 
asked to enler the starling and ending sample of the template in order 

15 to construct tho first template. Later, the template can be changed if 
there is a need for IL It can be done manually in the same way as the 
first template. It can be done also automatically in that way that the 
predetermined amount or lalesl delected P-waves is taken and 
averaged, wherein the averaged P-wave is set to be the next template. 

20 The lemplale is also checked for the adaptation, if the amount of 
dropped P-waves exceed the predetermined level. In these ways, the 
template can be changed according to the changes of the P-wave. In 
other words, template changes as the P-wave changes. Obviously, the 
P-wave can also be detected by some other methods known in the 

25 prior art, such as thresholding, patter recognition etc. 

It has been empirically observed that the preprocessed FCG-signal can 
still be relatively noisy around the P-wave, which can affect the 
determination of the real location of the onset and offset of P-wave 

00 The amplitude of the deflection might be affected by the existence of a 
noise plateaux before the beginning of the P-wave. For avoiding this an 
improved method of detection of the P-wave time limit is preferably 
used (309). in the improved method, the rising and falling edges of the 
P-wavc arc approximated by straight linos using a slope function. Later 

35 values of the baseline of the beat are calculated and stored 
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Intersection points of the baseline and the two edges are obtained by 
an intersection function. 

With the help of the baseline, the value of The onset and Die onset will 
5 be determined (310). The onset is the intersection of tho baseline and 
the rising edge, whereas the value of the offset will be the intersection 
of the baseline and the falling edge. The value for tho P-wave offset 
can also be found by a threshold method, in which the sample having 
the ECG of approximately tho samo value as in the P-wave onset is 
10 searched. The found values is accepted as the real offset of Ihe signal 
based on the idea that in most of the cases the P-wave loop will be 
closed in three-dimensional environment, which means (in lime 
domain) that the P-wave must return to its initial value. 

15 The method of finding tho P-wave onset and offset is extremely 
Important in analyzing the P-wave specific parameters, because in 
those calculations tho P-wavo must have clear boundaries, in other 
words, a P-wave vector loop must be completed. The P-wave vector 
loop is discussed in more detail later in the description. However, it is 

20 equally Important to notice that the offset function is not used while 
analyzing the PQ-eegment specific parameters, because with those 
parameters the PQ-segment elevation, in other words the incomplete 
P-wave vector loop, must be analyzed in the unedited situation. 
Therefore, the P-wave onset function is used in the method 

25 substantially all the time, and the H-wave offset function is used with 
tiie P-wave specific parameters. 

The detected P-wave is stored in X-, Y- and Z-leads (Hgure 4a, P- 
wave circled; Figure 3, 311). All the P-wave beats are then averaged 
30 (312) in the predetermined time intervals to form a smooth beat. 
Formed beats are used for calculations and analysis. The first 
averaged P-wave is used as an initial reference P-wave, where ihe 
upcoming averaged P-wave is compared ro (monitoring of the dynamic 
changes). 

35 
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When monitoring the myocardial infarct, tho patient is in cardiac 
monitor, e.g. MIDA, preferably for 1-2 days. Substantially all the ECG- 
signals during that time should be taken for the analysis. Averaging for 
the data Is preferably done every 4 minutes 

5 

When monitoring the cardioversion of the patient with atrial fibrillation, 
the normal sinus rhythm is observed prcforably for 3-4 hours 
Averaging for the data acquired Is preferably done every 10 sec - 2 
minutes. The reason for tho shorter averaging time is due to the lesser 
1 0 data and the quicker changes in the beginning. 

It is also possible to monitor the P-wave changes, for example. 1 3 or 
6 months after the procedure such as, for example, the cardioversion 
In that case Ihe monitoring may last 15-30 minutes and it focuses on 
15 observmg the change of the vector loop. The different monitoring 
results are combined for belter analysis. 

It Is obvious that the monitoring depends on ihe situation, when also 
^ the above-mentioned time limits vary according to the situation. 

In some situations, the P-wave can be seen as two loops, the primary 
' secondary (A„) loop, which are detected (313). The loops 

are Illustrated in figures 4b and 5b in one-dimensional and two- 
d.mensional spaces. Detection of the secondary loop Is essential in P- 
25 wave loop analysis, as it can be greatly valuable in the analysis of atrial 
anatomic changes described later in the text. P-wave loops in a three- 
d.mens.onal space can preferably be projected onto the three 
orthogonal planes and the existence of the secondary loop can thus be 
completed in two dimensions 

30 

The QRS-comp, ex is also observed in ^ metnQd accord|ng ff> ^ 
jnvention. Th.s is done by detecting the beginning and the ending of the 
QRS-complex from the beginning of the bear. The duration of the QRS- 
complex .s also measured, as well as the waveform signal. These 
35 values are used along the P-wave analysis 
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According to the invention, tho averaged P-waves are calculated and 
the parameters described later In the text are estimated next. (314). A 
few of the parameters arc common parameters, several are developed 
for the Invention and for the dynamic analysis of the P-wave. since they 
5 have not existed before. Parameters can bo divided into four 
categories: one-dimensional environment, two-dimensional 
environment, three-dimensional environment and magnitude ECG 
environment. 

10 One-dimen sional anvirnnm^ pt 

The parameters of a one-dimensional diagram describe Ihe P-wave 
properties in one dimension and relate to the orthogonal axes of the 
electrocardiogram. Figures 4a and 4b illustrate the P-wave in a one- 
15 dimensional environment A vector area (P-A) (Fig. 4a) expresses the 
area of Ihe P-wave. II can be derived from the areas in the three 
orthogonal leads X, Y. 2 to form an equation (sqrt corresponds to the 
square root): 

20 P-A = sqrt ( A/ + A Y 9 + A/ ). 

A vector change area (PC-A) describes the changes of the P-wave 
area. Preferably the area of the P-wave being examined (A^,,) is 
compared with the area of the reference P-wave (A, ftf ). The difference 
25 is calculated from X, Y, 2 -leads: 

PC-A = sq rt [(A^ - A^ H (A^ - Arcf ) Y ^(A^ - A fef ), \ 

A P-area duplicity (P-AD) can be calculated by relating the area of the 
30 secondary loop of the P-wave (A„) to the primary loop area (A,) from X, 
Y, Z-leads (Figure 4a). 

P-AD = sqrt[(A n / A,) K 2 + (An/ A|) Y a +(A| ( / A|) z 2 ]. 

35 In addition, the PQ-parameters are calculated. A PQ-vector magnitude 
(PQ-VM) describes the PQ-magnitude (level of PQ segment elevation) 
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and a PQ area (PQ-A) describes the area between the PQ-segment 
and the haseline (e) in a one-dlmensional environment (Figure 4c). A 
PQC-A-parameter describes correspondingly the upcoming changes of 
the PQ-area compared to the values of the reference P-wave. A PQ- 
5 time and a P-wave duration (P dur) are also calculated in the one- 
dlmenslonal environment. 

PQ-VM = sqrt( PQ X 2 +PQ Y 2 +PQ E 2 ) 
PQ-A = sqrt( PQA X 2 +PQA Y 2 + PQAz 2 ) 
10 PQC-A= sqrt [(A e x a m-A ref ) x 2 +(A < , xam -A ro ,) v 9 +(A exanl -A ref ) z 2 ] 



Two-dimensional environment: 

15 The parameters of a two-dimensional diagram describe a P-wave 
vector loop in the three orthogonal planes: frontal (XY), horizontal (XZ) 
and sagittal (YZ) planes. The P-wave vector loop in a two-dimensional 
environment is shown In figure 5a and 5b. A P-wave vector loop area 
(P-LA) is the area of a two-dimensional loop described by the P-wave 

20 vyelor (Figure 5a). II can be derived from the different areas in the 
three orthogonal planes XY, XZ and YZ to torm an equation: 

P-LA - sqrt ( A^+A^+A^ 2 ) 

25 A P-wave change loop area (PC-LA) is the difference of the two- 
dimensional areas of the P-wave vector loop being examined and the 
reference P-wave vector loop: 

PCLA= sqrt [(A 8Xam -A ref ) XY a +(A wam - A ref ) X2: 2 +(A exom -A rof ) YZ 2 l 

30 

A vector loop area duplicity (P-LAD) describes the P-wave duplicity 
ratio, which is the ratio of the secondary loop (A,,) to the primary loop 
(A) (Figure fib): 

35 P-LAD = sqrtr (A„ / A,) xy % <A„ / A,) xz 2 +(A„ / A,) yz ? ] 
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Three-dime nsional environment- 

I he parameters of a Three-dimensional diagram describe the P-wave 
vector loop in the three-dimensional environment (Figure 6a). A thrco- 
dimensionai P-wave vector loop area (P3-I.A) Is the area of the three- 
dimensional loop described by the P-wave vector. The basic function, 
called "looparea", calculates the area of a loop dividing It into small 
triangles in three dimensions and using the known vector cross product 
method to compute their areas. The three-dimensional P-change loop 
area (P3C-LA) is the difference between the area of tho P-wave vector 
loop being examined (A eKmn ) and the reference P-wave vector Iood 
(Am): 

15 P3C-LA = sqrt[(A eXam -A fet )% 

A P-azimut (P-Az) is a parameter for the angle that a P-main vector (M) 
describes In the transversal plane (Figure 6b), The P-main vector is 
defined as the average vector of substantially all the P-wave vectors 
PO that compose a P-wave loop in three dimensions. A P-Elevation (P-EI) 
shown also in Figure 6b, is a parameter for the angle that the P-main 
vector forms in the vertical plane. 

A P-change veclur angle (PC-VA) denoles the angle difference 
25 between tho main vector of the P-wave examined and the reference P- 
wave. A P-VM-pararneler is ihe P-veoLor magnitude and a P-change 
vector magnitude (PC-VM) is the magnitude of the vector point ot the 
reference P-main vecLor to the examined P-main vector. A P-QRS- 
«n X e ol 0r an9 ' 9 (PQRS ' VA > ^presses the angle between the P and the 
30 QRS-main vectors. A P-QRS-change vector angle (PQRSC-VA) 
expresses the upcoming changes compared to the values ot the 
reference P- and QRS-waves. 

t PQ " Vector ma » nftu * (PQ3-VM) describes the PQ-magnitude in a 
P^TT "^nnmwr. a PQ-chsnge vector magnitude 
(PQ3C-VM) .s the d.stance of the vector point of the reference P main 
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vector to the examined P-main vector magnitude in a three- 
dimensional environment. 

In addition, angles of the PO- vector (PQ-Az, PQ-FI) and the difference 
5 (PQC-VA) between the PQ-wave examined and the main vector of tho 
reference PQ-wave are calculated. 

A P-vector loop length (P-VLL) is the perimeter of the loop drawn by 
the P-wave vector in a three-dimensional space. A P-vector loop 
1 0 velocity (P-VLV) describes the speed of the development of the P-wave 
vector loop length. 

Magnitude ECG; 

15 The parameters of the magnitude ECG are a P vector magnitude area 
and the difference of it. Figure 7 illustrates the P-wave in magnitude 
ECG. The P-vector magnitude area (P MA) is an area of the P-wave in 
the magnitude ECG-signal. A P-vector magnitude area difference (PC- 
MA) describes the ohange between the examined P-wave magnitude 

20 area and the reference P-wave magnitude area: 

PC-MA=,sqrt[(A oxam -A rof ) ? ]. 

In addition, ihe PQ-paramelers are calculated similar to the one- 
25 dimensional environment A PQ-MVM describes the HQ-magnitude 
and a PQ-MA describes the PQ-area of the PQ-elevation. PQC-MA 
describes correspondingly the upcoming changes of the PQ-area 
compared to the values of the reference P-wave. 

30 PQC-MA = sqrt (A exam - A, e f) 2 

Figure 8a represents the change of the P-wave area (PC-A) and figure 
8b represents the change vector angle (FC-VA), as well as the change 
vector magnitude (PC-VM). 

35 
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Parameters described above are preferably also used with the analysis 
ot the beat stored Into the atrium extra systole database. 

When parameters have been calculated, the results are displayed 
(215). According to the invention, the results are displayed from every 
averaged time interval as a new point in trend curve over time. An 
advantage of this presentation is that the results arc easier to see and 
a conclusion is easier to draw from the changes of the P-wave. 

10 Analysis nf the parameter*' 

Parameters that concern the Pa-segment can have a great impaul 
when studying the development of an acute myocardial infarction of 
atria. It is known that the manifeslaliuns of atrial infarction may include 
15 elevation or depression of the Ta-segment (Ta represents atrial 
repDlarisation like the ST-segmeril and T-wave represent ventricular 
repolarisation). The Ta-sogmcnt is usually obscured by the qrs- 
complex and the early part of the ST-segmenl. but it's abnormalities 
may affect the PQ-segment. The parameters PQ-VM, PQ-A, PQC-A 
PQ3-VM, PQ3C-VM, PQ-MVM, PQ-MA. PQC-MA, PQ-Az, PQ-EI and 
PQC-VA are developed to analyze these changes. 



20 



25 



30 



35 



Parameters that concern the P-wave primary (A,) and secondary (A.,) 
loop are developed most of all to evaluate the P-wave changes in heart 
insufficiency, it is known that in heart insufficiency the P-wave may 
become peaked or notched, depending on the etiology of disease 
Although the specificity ot these changes varies, it is a common opinion 
that in the left atrial enlargement due to heart insufficiency the P-wave 
became notched and has a negative terminal part, so that the P-wave 
has the appearance of the "fallen S'\ The size of a negative terminal 
part may correlate to the level ot the left atrial enlargement and heart 
insufficiency. The parameters P-AD and P-LAD are developed to 
analyze these dynamic changes. 

Parameters that concern the QRS-complex vector angle and its 
ro!at.on to the P wave voctor angle arc developed to discover if the 
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change of P-wave vector angle is due to change of the heart position. 
The parameters PUHS-VA and PQRSC-VA are developed to analyze 
these changes. 

5 Most of the parameters that concern the P-wave and PQ-sogmcnt 
area, the vector magnitude and angle, the PQ-tlme and the P-duration. 
the length and the velocity of the vector loop describe the character of 
P-wave vector loop widely. Their changes reflect the condition of the P- 
wave dynamically. For example, in the acute myocardial infarction of 
10 atria these parameters will change due to the tissue damage of the 
atrium. 

There is some evidence that in the 12-ioad ECG the P-wave duration 
and dispersion (difference between the longest and shortest P-wave 

1 5 duration) may be greater in persons who arc in a higher risk to develop 
atrium fibrillation. Although the specificity of these changes varies, 
there is some evidence that when these changes decroase with the 
medication, the risk also decreases. As these changes in a 3- 
dimcnsional environmont can be seen as increased irregularity of the 

P0 P-wave vector loop, the above-mentioned parameters describe these 
dynamic changos ideally. Parameters are suitable also in other 
situations where dynamic changes in inlra-atrial pressure, volume or 
conduction alter the P-wave configuration slowly. In addition, these 
parameters are useful in investigation of the quality atrial extra 

25 systoles. 

System: 

A cardiac analysis system according to the invention comprises means 
30 for signal processing and counting parameters. The system is adapted 
to take raw ECG-signal from existent datasystem and to store it into a 
file. Said cardiac analysis system is arranged to convert the sample file 
into a binary file containing X, Y, Z -samples. The samples in the 
generated file are preferably 16 bits signed integers. Samples are 
35 advantageously ordered in X(i), Y(i), 7(i). X(i+1), Y(l+1), Z(l+I), wherein 
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X, Y, Z refers to the three orthogonal components of the VCG-signal 
and T refers to the beat number. 

An Input-file also .comprises Information about the resolution of the 
5 system, samples taken per second, data type (preferably signed 1 6 bit 
integer), beginning of the n me frame t0 be ^a^zea related lo Ihe 
whole duration of .the file, end of the time frame to bo analyzed and 
number of channels (preferably three, corresponding to the orthogonal 
axes X, Y, 2). Naturally, the input-file can include other information as 

1 0 Wfill. 
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The system also comprises means for preprocessing lo remove the 
no.so of tho signal and to improve the SNR. Preprocessing is 
preferably done by Tillering, which can be implemented with both low 
and h,gh frequency of the bandpass preprocessing filter. The system 
also comprises a structure, which focuses on the beat specifications. 
Preferably, specification is done by measuring time and duration ot tne 
peaks and the waves. II can be adapled to measure the previous R- 
peak the posterior R-peak, as well as the beginning ot the range to 
20 calculate ihe baseline. Also, il can be slored with information of the 
channel used for R-peak detector and baseline correction, the rate to 
calculate the threshold for R-peak detection, and the duration of every 
beat in samples. y 

Ihe detected beat is stored into the system. Information of the beat 
such as a samptenumber when the T:th beat starts in the entire input 
file; a sample number of the ending of the T:th beat in the input file- 
duration of the beat in samples; wave form signal of the stored beat 
basel.ne value of the T:th beat; duration of the beat in seconds;' 
location of the R-wave peak in samples; validity of the P-wave (0 or 1) 
deterrnmed preferably by duration; P- wave onset; P-wave offset and 

™2L P " WaVe Samp,eS are P referab 'y stored into the 

system. 



25 
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WaV l de T^ ° f the System fs capable ° r validating the P-wave by 
know.no the threshold of covariance of the template method and the 
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minimum and the maximum P-wave durations in samples. The validity 
of the P-wave Is determined by measuring P-wave duration, which 
should settle between tho minimum and the maximum values. 

5 The system is preferably implemented as hardware and software. 
Therefore the system also comprises a computer program for 
implementing the ^method according to the invention. The computer 
program comprises computer readable code for acquiring the ECG- 
data, detecting the P-wave from said data and analyzing said P-wave 
10 according to the method of the invention. 

The method and the system describe the preferred embodiments of the 
P-wave analysis according to the invention. The main idea in the 
method is to analyze the dynamic changes of the P-wave in time. 
15 Implementation of :the system can be carried out in different ways. By 
knowing this, it is obvious that the present invention is not limited to the 
description, but to claims discussed herein below. 
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1 . A method for a cardiac analysis, the method comprising steps for 
acquiring an ECG-siynal, deLecling at least one wave of the ECG- 
signal and calculating parameter values of said wave, characterized in 
that, said wave is a P-wave, whereupon the cardiac analysis is focused 
to dynamic changes of the P-wave, wherein substantially every 
detected P-wave is compared lo a reference P-wave in defined time 
period. 

2. The method of ; claim 1, characterized in that the bCG-signal is in 
the form of a vectorcardiogram. 

3. The method of :claim 1 or 2, characterized in that a beat between 
two R-peaks is examined, whereupon said beat is classified into 
groups depending on whether the beat is having a duration between 
the predetermined time limit or the beat is having a duration under the 
predetermined time limit, whereupon both said beats are analyzed 
separately. 



claims 1 - 3, characterized in that the P-wave is 



35 



4. The method of 
detected by a template method. 

5. The method of- claims 1 - 3, characterized in that the P-wave is 
detected by a pattern recognition method. 

6. The method of ; clalms 1 - 5, characterized in that the detected P- 
wave is stored in X, Y, 2 Iead3. 

7. The method exclaims 1 - 6, characterized in that the detected P- 
wave Is averaged jin the predetermined time Interval. 

8. The method of plaims 7, characterized in that at least one averaged 
P-wave is used as an initial reference P-wave, where the upcoming 
averaged P-waves are compared to. 
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9. The method of claims 1 - 8, characterized in that at least one loop 
of the P-wave is detected. 

10. The method o{ claims 1 - 9, characterized in that the parameters 
5 of the p-wave lnj one-dimensional diagram are one or more of the 

following: the vector area (P-A), vector change area (PC-A), P-area 
duplicity (P-AD), PU-vector magnitude (PQ-VM), PQ-area (PQ-A) and 
PQ change area (PQC-A). 

10 11. The method ofjclaims 1-10, characterized in that the parameters 
of the P-wave in i two-dimensional diagram are one or more of the 
following: the vector loop area (P-LA), vector change loop area (PC- 
LA) and P loop area duplicity (P-i AD). 

15 12 The method of: claims 1 - n . characterized in that the parameters 
of the P wave in three-dimensional diagram are one or more of the 
following: the vector loop area (P3-LA), the vector change loop area 
(P3C-LA), the angjos of the azimun (P-Az, PQ-Az), the elevation (P-EI 
PQ-E1). change vector (PC-VA, PQRSC-VA, PQC-VA) the P-QRS 
vector (PQRS VA); as well as the vector magnitude (P-VM, PQ3-VM) 
change vector magnitude (PC-VM, PQC3-VM). 

1 f'J hB o eth ° d of : :claims 1 ~ )2 > characterized In that the parameters 
ot the P-wave in magnitude environmont arc one or more of the 
following: the vector magnitude area (P-MA, PO-MA). the vector 
change magnitude area difference (PC-MA, PQC-MA) and tho vector 
magnitude (PQ-MVM). 
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14. The method df claims 1 - 13, characterized in thai Lhe method 
compnses also the; calculations of one or more of the following: the PQ- 
time, P-wave duration (P-dur), the ler.ylh of the P-wave (P-VLL) the 
velocity ot the P-wkv© vector loop (P-VLV). 



35 
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15. The melhud of claim 1 - 14, characterized in that the ECG-signal 
ie acquired from a Frank system or a 12-lead ECG-arrangement 

16. The method exclaim 1 - 14, characterized in that the ECG-signal 
is acquired from a .JWDA data storage unit. 

17. The method of claims 1-16, characterized in that results of the 
parameters are displayed In a trend curve. 

18. A cardiac analysis system comprising tirst means for acquiring the 
ECG-signal, second means for detecting at least one wave from the 
ECG-signal and third means tor calculating parameter values of said 
wave, characterized in that, said wave is a P-wave, whereupon the 
cardiac analysis system is adapted to tocus to dynamic changes of the 
P-wave, wherein r said system additionally comprises means for 
comparing substantially every detected P-wave to a reference P-wave 
in defined time period. 

19. The system of;claim 10, characterized in that the ECG-signal is in 
form of a vectorcardiogram. 

20. The system of claim 18 or 19, characterized in that the system is 
also adapted to measure a duration of the beat between two R-peaks 
wherein the system is also configured to compare the heat to the 
predetermined time limit and classified the beat into the one of two 
groups depending on whether the duration Is between the 
predetermined time limit or under the predetermined time limit, wherein 
the system is alsoponfigured to analyze both groups separately. 

21. The system ofi claims 18-20, characterized in that the system is 
adapted to detect tho P-wave by a template method. 

22. The system exclaims 18 - 21, characterized in that the system is 
adapted to detect the P-wave by a pattern recognition method. 



VAST.0TT0 1D-Dd-ZDQ3 16;07 MISTA- 03 ZB8BZBZ 



KENELLEPATREK Asiakaspalval SIVU 023 



23 



10 



15 



20 



25 



30 



35 



23. The system declaims 18-22, characterized in that the system is 
adapted to store trie detected P-wave in X. Y, Z leads 

24. The system exclaims 18-23, characterized in that the system is 
adapted to average the detected P-wave in the predetermined time 
interval. 

25. The system offclaims 18-24, characterized in that the system is 
adapted to use the. first averaged P-wave as a reference P-wave and to 
compare the upcoming averaged P-waves to it. 

26. The system of?claims 18-25, characterized in that the system is 
adapted to detect at least one loop of the P-wave. 

27. The system ofjclaims 18-26, characterized in that, the system is 
adapted to acquire the ECG-data from a Frank system or a 12-lead 
ECG-arrangementi 

28. The system ofjclaims 18-27, characterized in that the system is 
adapted to acquire the ECG-signal from a MIDA register. 

29. The system exclaims 18-28, characterized in that the system is 
adapted to displayjresults of the parameters calculated in trend curve. 

30. A computer program product, comprising a computer readable 
storage medium on which is stored a computer program code for a 
cardiac analysis,;- which computer program code comprises first 
computer instructions configured to acquire tho ECG-signal, second 
computer Instructions configured lo delect at least one wavy from the 
EOb-signal and third computer instructions configured to calculate 
parameter values ;pf said wave, characterized in that, said wave is P- 
wave, whereupon; the computer program code has instructions for 
focusing to the dynamic changes of said P-wave, wherein said 
computer program code additionally comprises computer instructions 
configured to compare substantially every detected P-wave to a 
reference? P-wave In defined time period. 
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(57) Abstract: (_3 

The invention relates lo a method, lo a system and lo a 
computer program; product for the cardiac analysis. The 
method comprises|steps for acquiring the ECG-signal in 
the form of vectorcardiography, detecting P-wave from 
the ECG-signal by a template method and calculating 
the parameter values of the P-wave preferably 
continuously during the ECQ-recording. The method is 
aimed to the dynamic changes of the P-wave, wherein 
substantially every: detected P-wave is compared to the 
reference P-wave \ in defined time period. I he cardiac 
analysis system according to the invention is configured 
to implement the aforementioned method. I he computer 
program product comprises computer program code for 
implementing the aforementioned method. 

(Fig. 3) \ 
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